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SUMMARY 

/ 6^ 6 

Detonation induction distances were obtained experimentally 
for hydrogen- oxygen mixtures having fuel concentrations of U 5 , 
66.67 and 75 per cent hydrogen at initial pressures of 1 and 5 
atmospheres in a detonation tube with variable internal geometry. 
When three stainless steel rods of different diameters were 
located concentrically in the initiation region of the tube, the 
induction distances decreased as the diameter of the rods was 
increased. With a divergent, stepped-wall insert placed in the 
initiation section of the tube, the induction distances were 
reduced greatly compared to those obtained in a straight cylin- 
drical tube at 1 atmosphere initial pressure. When the tube had 
a right-angle section located 66 cm from the ignitor, a detona- 
tion wave formed immediately past the corner for all fuel con- 
centrations at 1 and 5 atmospheres initial pressure except for 
the stoichiometric mixture at 5 atmospheres pressure for which 
a detonation wave formed in less than half the distance from the 
ignitor to the corner. A detonation wave formed again at or 
prior to the right-angle corner, when it was located 15 cm from 
the ignitor, at initial pressures of 1 and 5 atmospheres except 
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for ^5 and 75 per cent mixtures at 1 atmosphere pressure where 
the induction distances were greater than those derived from 
the other position of the corner; however, these distances were 
shorter than those obtained in a cylindrical tube. Increased 
turbulence and the generation of stronger pressure waves from 
the various geometrical configurations are believed to cause 
the significant decrease in detonation induction distances in 
variable geometry tubes of the type investigated. 


When a combustion wave reaches an area discontinuity or 
a rod-insert, pressure waves will be generated. These second- 
ary pressure waves will lead to wave-interactions which will 
strongly influence the transition process from deflagration 
to detonation. On the basis of the small perturbation analysis, 
it is shown that the effects of the area disturbances generally 
will enhance the transition process, as verified by measurements 
of the detonation induction distances. 


The results of the measurements of detonation induction 
distances with various rod- inserts and geometrical configura- 
tions are included in this analysis. It is evident that the 
mechanism of turbulence generation is quite important when 
these obstacles are placed into the flow field. Without con- 
sideration of the effects of turbulence, the transition mech- 
anism for a f lam e propagating in a duct with nonuniform cross 
section cannot be explained satisfactorily. 


EXPERIMENTAL INVESTIGATION 


A. Introduction 

In the study of combustion phenomena, the transition of a flame 
from deflagration to detonation arises as one of the more important 
problems. Complicating the study of this process is the fact that 
the transition mechanism is an unsteady process involving flame 
acceleration, shock-deflagration-rarefaction wave interactions, gener- 
ation of turbulence and effects resulting from chemical and molecular 
kinetics. The distance and time required for the flame to accelerate 
to the detonation state is long for some mixtures and very short for 
others. Many factors influence the distance which is required to 
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establish the detonation wave from the initial deflagration flame. 

Some of the parameters known to have a noticeable influence on the 
formation of the detonation wave are the initial pressure, initial 
temperature, fuel concentration, type of ignition, linear flow 
velocity of the unbumed gas in a tube, the variation in internal 
geometry of the tube containing the gas mixture and external magnetic 
fields (Refs. 1-7). 

A practical parameter describing the rate of formation of a 
detonation vave is the detonation induction distance. As used in this 
report the detonation induction distance is the distance between the 
source of ignition and that point in the detonation tube where the 
propagation rate first attains a velocity which is equal to the stable 
detonation velocity (Chapman-Jouguet state). Interest in the formation 
of detonation waves arises from the technical and fundamental aspects 
of combustion. Results from the studies of the formation of detonation 
waves are expected to contribute to the understanding of combustion 
instability in rocket engines and to the development of engines employ- 
ing detonative combustion. 

The experimental portion of this investigation was conducted to 
determine what effect different geometric configurations of the deto- 
nation tube would have upon the detonation induction distances for 
various hydrogen-oxygen mixtures. With each tube configuration, 
numerous experiments were made at different fuel concentrations and 
initial pressures. Various obstacles were placed into the tube. To 
induce turbulence in the center of the tube and to change the cross- 
sectional area of the region of flow passage, three stainless steel 
rods were inserted concentrically in the initiation section of the 
detonation tube. In effect, the arrangement created an annular tube 
in this region. The turbulence in the expanding gas was increased 
by these inserts . 

To study the effect of a sudden expansion of the tube cross- 
sectional area upon the induction distance, a divergent, stepped-wall 
insert was employed. With this arrangement the gas expanded suddenly 
when it reached the step in the insert. While the linear gas 
velocity decreased suddenly at these steps which reduced the formation 
of a detonation wave, the amount of turbulence was increased due to 
flow separation and eddies at the steps. 

A 90-degree angle in the initiation section of the detonation 
tube permitted a study to be made of the effect that angular motion 
of the gas has upon the establishment of a stable detonation wave . 

In one arrangement the 90-degree corner was located only 15 cm from 
the ignitor while in the other set-up the corner was 66 cm distance 
from the source of ignition. 


0 
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The results of all these experiments give an indication of the 
effects of these three major variations in tube configuration. 

Although there are many other possible geometric variations which 
should be investigated to permit a thorough understanding of the 
influence of tube geometry on the formation of detonation waves, it is 
believed that some of the most important ones have been selected and 
investigated. 


B. Experimental Equipment 

All of the experiments concerning the formation of detonation 
waves were conducted in a stainless steel tube with a 79-mm inside 
diameter and a wall thickness of about 2 cm; a photograph of this 
tube is shown in Fig. 1. For the measurements in a straight cylin- 
drical tube, four one -meter long flanged sections were bolted 
together. An exploded view of one section is depicted in Fig. 2. 
Electric heaters, having 23 kilowatts capacity, were mounted exter- 
nally on the tube sections and flanges in order to maintain the tube 
at a selected initial temperature. Thermocouples and an electronic 
temperature regulator system maintained the detonation tube at a 
fixed temperature of 40°C for all experiments. To maintain this 
temperature as accurately as possible, one -quarter-inch diameter 
copper tubing was wrapped around each section of the tube so that 
water flowing through the tubing provided an additional heat sink 
thereby improving the thermal lag characteristics of the control 
system. 

The fuel and oxidizer gases were taken from standard cylinders 
and passed through a manifold filter to a silica gel trap to remove 
any oil or water vapor (Fig. 3)* The silica gel trap was maintained 
at a low temperature by the use of a freon refrigeration system. 

After metering, the two gases entered the mixing vessel which was 
located close to the detonation tube both of which were set up in 
a firing pit. Appropriate valving was used so that the detonation 
tube could be purged with dry air, evacuated and then filled with 
the explosive mixture of gases. When initial pressures higher than 
atmospheric were needed, the tube was vented to the atmosphere through 
a suitably adjusted back-pressure regulator. This procedure of mixing 
the gases as they were needed avoided the problems and hazards of 
filling a high-pressure supply tank with a large amount of uniformly 
premixed gas for a series of experiments. Valves which were employed 
to seal off the detonation tube and to vent the premixed gas were 
operated remotely because of the danger of a possible explosion. 
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Probes utilizing the conduction properties of the wave were used 
to detect the passage of the combustion or detonation wave as illus- 
trated in Fig. 4. A Teflon -coated wire was inserted through the wall 
of the tube so that the probe wire was flush with the inside surface 
of the tube. A 22. 5 -volt battery and a coaxial line to a multi- 
channel chronograph were connected in series between the wire and 
ground. When the combustion or detonation wave passed the probe, the 
ionic content of the wave was sufficient to trigger the chronograph. 

The time intervals between various probe positions were determined 
by the difference in the chronograph readings to the respective probes . 
A six-channel ten-megacycle chronograph and a ten-channel one -mega- 
cycle chronograph, Fig. 5 > which were developed in this laboratory, 
were utilized for these measurements. Details of the chronograph 
system can be found in Ref. 8. Probe positions were spaced along the 
tube at numerous places to provide adequate coverage for the velocity 
measurements . The type of ignitor employed for all of the experiments 
consisted of a 3/8 -inch length of 0.005 -inch diameter copper wire 
which was melted with current from a 28-volt D.C. power supply. 


C. Tube Geometry 

In the study of the effect of tube geometry on the detonation 
induction distances, three cylindrical rod inserts of different diam- 
eters (Figs. 6-8) and a diverging insert (Fig. 9) were located near 
the ignitor end of the detonation tube. The three stainless steel 
cylindrical bodies were mounted concentrically within the detonation 
tube; they extended approximately 84 cm downstream from the ignitor. 
The diameters of the bodies were 0.635 cm , 1.59 cm and- 4.76 cm, and 
the ratio of the diameter of the bodies of revolution to the deto- 
nation tube diameter were 0.08, 0.20 and 0.60 respectively. 

The insert containing the divergent stepped-walls (Fig. 9) was 
machined from aluminum stock and had an overall length of 57*42 cm. 
One end of the insert was located at the ignitor station. This 
insert provided an increase in diameter by a factor of approximately 
6 in two abrupt steps. The first section of the insert, having a 
diameter of 1.35 cm, extended 24.13 cm downstream from the ignitor. 

At the end of the first section, the diameter increased to 3-81 cm 
which was constant for the rest of the length of the insert. The 
diameter then increased to 7*9 cm at the end of the insert which 
was the diameter of the detonation tube. Holes were drilled for 
the detection probe wires after the insert was assembled in the 
detonation tube . 
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An additional section which was fabricated for the detonation tube 
is shown in Fig. 10. It was constructed with a right -angle bend in the 
tube in order to study the effect of angular motion of the gas on the 
formation of detonation waves. The section was made from the same 
stainless steel stock that was used to fabricate the other four one- 
meter long sections of the detonation tube. 


D. Experimental Procedure 

In general, the same procedure was used for all experiments . 

After the heaters on the detonation tube were turned on, and the 
temperature of the tube was raised to 40°C, it was evacuated by the 
use of two conventional vacuum pumps. Evacuating the detonation tube 
shortened the flow time needed to assure a homogeneous mixture of 
hydrogen and oxygen in the tube, and the procedure also minimized the 
amount of contaminant gas. To fill the tube, the control valves on the 
flow stand were adjusted so that the desired flow rate for the partic- 
ular fuel concentration required was obtained. Next, the premixed gas 
from the mixing chamber was allowed to flow into the detonation tube 
until the desired pressure level was obtained by employing a back- 
pressure regulator which allowed the gas mixture to flow out the 
vent when the pressure was above a predetermined level. After the 
explosive mixture reached the desired pressure in the tube, the gases 
were allowed to flow for approximately two complete changes of free 
volume in the tube to ensure a uniform mixture. Then, the vent and 
fill valves of the detonation tube were Closed and the flow of gases 
was terminated. After the gas had reached temperature equilibrium, 
excess gas was bled off until the required level of initial pressure 
was reached. During the time needed to reach temperature equilib- 
rium, the chronograph system was checked with the internal test 
system (see Ref. l). All of the electronic equipment including the 
electronic gates, oscillator, pulse counters, and so forth were 
tested before each firing of the detonation tube. Prior to igniting 
the mixture, the electronic counters were reset to zero. After the 
explosions, the tube was vented and flushed with dry air to remove 
the combustion products and the next experiment was started. 

During previous investigations, it was found that identical 
initial conditions did not give identical results in the initiation 
region of the detonation tube. It was established that these differ- 
ences were not caused by the chronograph system. Time measurements 
from the region where stable detonation occurred were quite repro- 
ducible. Therefore, in order to obtain usable results it was necessary 
to make a number of measurements at each location of the probes for 
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each set of initial conditions. The velocity data from each station 
were averaged and graphed as a function of distance from the ignitor. 
In addition, the maximum and minimum values of measured velocity were 
indicated too. The detonation induction distance, based on the 
maximum propagation rates, was determined from the graphs as the point 
where the curve through the values of the maximum propagation rates 
reached the theoretical value of the detonation velocity for the 
particular set of initial conditions. The induction distance based 
on the average velocities was determined similarly. In a sense the 
induction distances determined from the maximum propagation rates 
give a "safe" induction distance if one is concerned with the destruc- 
tive aspects of a detonation wave. Of course, the variation in 
propagation rates in the ignition region are probably not true 
velocity changes. The indicated variations could result from 
spinning flame fronts, tilted flame fronts and turbulent fronts 
which could have spikes of flame extending forward along the wall 
of the tube where the turbulence is generated by friction. The 
forward-protruding spikes of flame would trigger the detection 
probes; repeated firings with the same initial conditions would not 
give the same results because the formation of the flame spikes 
undoubtedly is a random process. 


E. Results 


Straight Tube Without Obstructions 

Some experiments were conducted with the straight, cylindrical 
tube before any obstructions were placed inside. Three fuel concen- 
trations, 45* 66.67 a^d 75 per cent, were employed at initial pressures 
of 1 and 5 atmospheres with hydrogen and oxygen. Graphs of the 
velocity variations are illustrated in Figs. 11 through l6. The 
induction distances are graphed in Figs. 27 and 28; the data are 
compiled in Table 1 along with the induction distances obtained from 
experiments with the cylindrical rod inserts. 

Since the diameter affects the induction distance, a comparison 
can be made between the data obtained with the 79 *™"* tube and data 
obtained during previous experiments with 15 -and 50 -mm diameter tubes 
(Refs. 1 and 2) as shown in Table 2. The available induction distances 
are listed for both maximum and average flame propagation rates at 1 
and 5 atmospheres initial pressure. No data were taken at 5 atmos- 
pheres initial pressure in the 50 -rnm diameter tube. 

From Table 2 it can be seen that the detonation induction distances 
become longer as the tube diameter increases for 45 and 75 per cent fuel 
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concentration at both 1 and 5 atmospheres initial pressure. At stoi- 
chiometric proportions, however, the induction distances are practi- 
cally the same in the 15- and 50-mm diameter tube; however, they are 
longer in a 79-mm diameter tube. There is no simple relationship 
evident from these data between detonation induction distance and tube 
diameter. Generally, the effect of tube diameter on detonation 
velocity can be expressed in terms of a linear relationship between 
detonation velocity and the reciprocal of the tube diameter. No such 
expression exists for the induction distances as the tube diameter 
changes. 

From these limited results, it does not appear that an empirical 
relationship can be developed to express the dependence of detonation 
induction distance on tube diameter for hydrogen and oxygen mixtures. 
More measurements are needed to elucidate this effect. 

Cylindrical Rod Inserts 

A description of the three rods used is given in Section C; a 
cross-sectional view of their location in the tube is shown in Figs. 

6 through 8. Experiments were conducted at fuel concentrations of 
45, 66.67 and 75 per cent hydrogen at initial pressures of 1 and 5 
atmospheres with an initial temperature of 40°C. All mixtures were 
ignited with the exploding wire ignitor. 

Experiments with the O.635- and 1. 588 -cm diameter inserts were 
made only with a fuel concentration of 66.67 per cent hydrogen at 
initial pressures of 1 and 5 atmospheres. Graphs of the average, 
maximum and minimum flame propagation rates as a function of distance 
from the ignitor are shown in Figs. 17 through 2 6. Also shown on 
these graphs are the induction distances as determined from the curves 
of the m aximum velocities and the average velocities. These induction 
distances are dependent upon fuel concentration as shown in Figs. 27 
and 28 respectively. The data are listed in Table 1 from which it can 
be seen that there is a definite decrease in the induction distances 
for both the O.635- and 1.588-cm inserts compared to the induction 
distances without an insert at 1 atmosphere initial pressure. At 5 
atmospheres, however, the induction distances for the 0.635“ and 
1.588-cm diameter rods decrease but not in the same proportion as 
those for the straight tube without obstructions which can be seen 
in both Figs. 27 and 28. The induction distances for these two inserts 
are practically equal to or greater than those for the straight tube 
without obstructions at 5 atmospheres even though the induction 
distances for the rod inserts are appreciable less than those for the 
same tube at 1 atmosphere pressure. 
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Thus, it can be seen that the effect of pressure on the induction 
distances is rather small when these rods were inserted in the tube as 
compared to the effect in a straight tube without obstructions . These 
results can be explained by considering the effect of turbulence which 
is introduced by the insertion of the rods. It is well known that 
increased turbulence enhances transition to detonation. It is also 
known that increased pressure reduces the induction distances of hydro- 
gen-oxygen mixtures in straight tubes without obstructions. Therefore, 
if the transition is forced to occur in a shorter distance by turbu- 
lence (e.g., produced by cylindrical rods), then it is to be expected 
that the effect of pressure on the transition distance will be reduced 
when compared to the results obtained with straight tubes without 
obstructions . 

During the experiments with the 0.635 -cm insert, very high veloc- 
ities (e.g., average velocities in excess of 10,000 meters /second) 
were obtained at a distance of 89.6-cm from the ignitor. These flame 
propagation rates were fictitious; therefore, they were not plotted 
on the graphs. The high apparent velocities were caused because of 
interference by the rod insert which was located 83.2-cm from the 
ignitor. If the flame front is non-planar, the detection probes on 
the wall of the detonation tube generate triggering signals to the 
chronograph which cannot be interpreted in the ordinary manner. Turbu- 
lence and shock waves from the upstream end of the rod disturbed the 
planar flame front. Probe signals could have been generated which 
gated the chronograph thereby giving rise to fictitious velocities 
because of the erroneous time interval. 

For the 4.763-cm diameter insert, the induction distances 
decreased greatly compared to the results of the experiments with the 
O.635- and 1.588-cm inserts. Again, similar to results obtained with 
the 0.635 -cm insert, there were some fictitious velocities obtained 
at different locations downstream from the ignitor. For an initial 
pressure of 1 atmosphere and a fuel concentration of 66.67 per cent 
hydrogen, fictitious velocities were obtained at two points located 
7.6 cm and 12.7 cm from the ignitor. Turbulence created by the forward 
end of the insert located 4 cm from the ignitor or the insert supports 
located at 9*08 cm from the ignitor could have caused erroneous time 
measurements. The same situation existed at 76. 1 cm for an initial 
pressure of 1 atmosphere and 75 per cent fuel concentration. Also, 
at this initial pressure and fuel concentration, a detonation wave 
was formed prior to the first insert supports (9.08 cm from the 
ignitor) for some experiments, but not in others. Evidently the 
amount of turbulence created by the probe mount on the upstream end 
of the insert was sufficient in some cases to form a detonation wave 
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earlier than in other experiments. Erroneous velocities also were 
obtained at 7*6 cm and 12.7 cm from the ignitor for a 45 per cent 
fuel concentration at 5 atmospheres initial pressure and at 12.7 
cm from the ignitor for a 75 P er cent fuel concentration at the 
same initial pressure. 


Divergent Insert 

In order to extend the study of how internal geometry affects 
the induction distances in hydro gen -oxygen mixtures, an insert 
(Fig. 9 ) containing divergent stepped-walls was employed. Experiments 
were made at an initial pressure of 1 atmosphere only; they were 
discontinued after the insert was damaged when the first mixture at 
5 atmospheres was detonated. The velocities obtained in the region 
of the insert varied considerably from experiment to experiment as 
shown in Figs. 29 through 31 * With 45 and 66.67 per cent fuel 
concentrations, the velocities obtained indicated that a detonation 
wave was formed in the first section (1.35 cm diameter) of the insert 
whereas at 75 P er cent fuel concentration no detonation wave was 
formed until the propagating flame reached the next section of the 
insert or was completely through the insert. Induction distances 
are listed in Table 3 * 

The effect of the insert on propagation velocity was quite evident 
from Figs. 29 through 31 * As expected, the velocity decreased in all 
cases when the gas expanded at each of the two stepped-wall sections. 
When the gas expanded to the tube diameter, 7*9 cm, the detonation wave 
collapsed for the 45 and 66.67 per cent mixtures. But, the flame 
propagation rates indicated that the detonation started to re-form, 
especially with the stoichiometric mixture (Fig. 30 ). The detonation 
wave did not collapse at the end of the insert with the 75 per cent 
fuel mixture. Only a modest decline in velocity was obtained (Fig. 31 )* 

Detonation induction distances determined from the maximum and 
average flame propagation rates in the divergent stepped-wall insert 
are compared with induction distances obtained in a straight tube 
without obstructions in Figs. 32 and 33 respectively. It is obvious 
from these illustrations that the induction distances with the insert 
are much lower than those obtained with the straight tube without 
obstructions. Since the 45 and 66.67 per cent mixtures detonated within 
the 1.35 -cm diameter portion of the insert, the induction distances 
should be comparable to those obtained with a cylindrical tube of this 
same diameter without an insert. Previously, data had been obtained with 
hydrogen-oxygen mixtures in a 1.5-cm diameter tube. With the insert, 
however, the detonation induction distances are lower by factors of 
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from 3 to 10 than those in the 1.5 -cm diameter straight tube without 
obstructions. The small difference in diameter from 1.5 to 1.35 cm 
does not appear to justify the drastic reduction in induction distance. 
It is believed that the large difference in induction distances results 
from the difference in vail roughness of the tubes. The 1.5-cm 
diameter tube had a very smooth vail vhereas the divergent steppea- 
vall insert was rather rough. It vas machined from aluminum, but it 
could not be honed. Therefore, it is believed that the resultant 
turbulence from the rough vail of the insert caused the detonation 
wave to form in a shorter distance. 

Right Angle Section 

To determine detonation induction distances with a right -angle 
section in the detonation tube, some experiments were conducted with 
the 90°-bend positioned 15 cm from the ignitor and additional exper- 
iments with the bend located 66 cm distant from the Ignitor. Initial 
pressures of 1, 5 and. 10 atmospheres were employed with fuel concen- 
trations of 45, 66.67 and 75 per cent hydrogen. The initial temper- 
ature vas 40°C. The propagation velocities obtained during these 
experiments are graphed in Figs. 34 through 51> and the induction 
distances are given in Table 4. 

At an initial pressure of 1 atmosphere with the 90° -bend 66 cm 
from the ignitor, a detonation wave vas formed Immediately after the 
corner for all fuel concentrations. The induction distances for these 
cases were estimated to be the same for both methods of determining 
the induction distances (maximum and average flame propagation rates). 
The same results were obtained with an initial pressure of 5 atmos- 
pheres with fuel concentrations of 45 and 75 per cent. At the stoichio- 
metric mixture, the induction distances (10 and 29 cm) were less than 
the distance from the ignitor to the comer. The induction distances 
in a straight tube for these Initial conditions are 37 and 56 cm based 
on maximum and average propagation rates respectively (Table 1). The 
exact reason for this discrepancy is not known. One possible expla- 
nation is that some additional turbulence was present in this partic- 
ular tube section. As anticipated from previous experiments, a deto- 
nation wave formed prior to the location of the corner for all three 
fuel concentrations when the initial pressure was 10 atmospheres. 

With the 90°-bend located 15 cm from the ignitor, the detonation 
wave did not form immediately at the corner for the 45 and 75 per cent 
fuel concentrations at 1 atmosphere initial pressure. At the stoichio- 
metric mixture, however, detonation occurred just downstream of the 
corner. Although turbulence was present at the corner, it appears 
that the amount vas insufficient to cause transition for the lean and 



12 


fuel-rich mixtures. Ordinarily, in straight tubes without obstructions, 
these mixtures require a length of the order of 200 cm to detonate. At 
15 cm distance, the flame propagation rates of these mixtures was not 
too fast and, therefore, the comer did not cause too much turbulence. 

When the initial pressure was increased to 5 and 10 atmospheres, 
the induction distances obtained were in the general region of the 
right-angle comer (15 cm from the ignitor). Only at the stoichiometric 
mixtures did detonation occur before the corner (5 cm using the maximum 
propagation rates). Since no data were taken with the straight tube 
without obstructions at 10 atmospheres initial pressure, no comparison 
can be made. It is reasonable, however, to assume that the induction 
distance at 10 atmospheres initial pressure is near 5 cm because it is 
known from other experiments with hydro gen -oxygen mixtures in smaller 
diameter tubes that the induction distance decreases rapidly with 
initial pressure for stoichiometric mixtures. 


F. Discussion 

Cylindrical rod inserts reduce the detonation induction distance 
in proportion to the diameter of the insert. The primary cause of the 
decrease in the induction distance is the Increase in the amount of 
turbulence which is created by the insert and its supports. With the 
divergent stepped-wall insert, the induction distance was reduced 
greatly in all experiments compared to the induction distances 
obtained for a straight tube without obstructions. Turbulence created 
by the sudden changes in area caused a reduction of the detonation 
induction distances . 

The induction distances obtained at 1 and 5 atmospheres initial 
pressure with the 90 ° -bend at 15 cm and 66 cm from the ignitor were 
less than those obtained with the straight tube without obstructions. 
For some initial conditions, the detonation wave formed immediately 
past the corner while for stoichiometric mixtures it did not. At 
the lower pressure, apparently the turbulence was enough to lead to 
detonation in a shorter distance than for the straight tube. The 
detonation wave formed prior to the 66 cm distance at higher initial 
pressures which was to be expected from previous experiments. With 
the right-angle section 15 cm from the ignitor and the mixture at 1 
atmosphere initial pressure, the distance that it took for the deto- 
nation wave to form varied considerably; these variations appear to 
depend on the fuel concentration and the amount of turbulence created 
by the corner. 
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Generally the overall trend with the various tube configurations 
employed in this investigation is a reduction in the induction distances 
compared to those obtained for a straight tube without obstructions. 


THEORETICAL INVESTIGATION 


A. Introduction 

The mechanism of transition from deflagration to detonation 
involves the complex phenomena of wave interaction, turbulence gener- 
ation, and chemical kinetics. The rationalization of any of these 
effects as being the predominant factor in the development of deto- 
nation waves provides a limited basis for theoretical analysis of 
this mechanism (Ref. 9). A comprehensive mathematical description of 
the inter-related phenomena leading to flame acceleration, and hence 
to the eventual formation of detonation waves, has not yet been explored 
successfully. In a different aspect, Bollinger, Fong and Edse (Ref. 

6) have studied the detonation induction distances both experimentally 
and theoretically for various combustible mixtures in a constant-area 
pipe. Despite the limited scope of this approach, the transition 
process in terms of the induction distance can be correlated empir- 
ically to an almost universal function involving the various physical 
causes for flame acceleration. By "almost universal" we mean that 
this function exhibits an almost unique relationship with the measured 
induction distances, at least for the particular constant-area pipe 
selected for the measurements . It has been pointed out also in Ref. 6 
that the empirical curve agrees fairly well with the experimental data 
for non-hydrocarbon systems but not for hydrocarbon systems . A tenta- 
tive explanation of this peculiarity is that the possible excessive 
ion concentration in the unburned medium ahead of the flame front for 
a hydrocarbon mixture may have actuated the ionization probes prema- 
turely, thus rendering the determination of the true induction distance 
much more difficult with this detection technique. 

If an area change or an obstruction is present in the combustible 
fluid medium in which a flame propagates, additional interaction effects 
will arise and will be superimposed on the already existing flame 
acceleration process. The pre compress ion wave will be partly trans- 
mitted and partly reflected upon reaching the area disturbance or 
impinging on the obstruction. The reflected wave will interact with 
the ensuing flame front while the transmitted wave will modify the 
unburned fluid medium. Thus, the flame will undergo a change in 
burning rate as well as in conf iguration, and will create compression 
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(or rarefaction) waves propagating both upstream and downstream of the 
flame front. As the flame front arrives at the cross-sectional area 
change or the obstruction, it will experience further changes in 
properties since a new set of boundary conditions for the flame front 
now prevails. Pressure waves emanating from the flame front will be 
generated whose type and strength are dictated by the particular 
disturbance the flame encounters. 

Besides the mechanism of wave interactions, the process of 
turbulence generation also is expected to be intensified under these 
circumstances inasmuch as distortion of the flame surface, vorticity 
generation, flow separation, and viscous drag on the obstruction, 
individually or collectively, will yield a strong influence on 
inducing flow and flame turbulence* Without considering the effects 
of turbulence, erroneous conclusions may be drawn sometimes. For 
instance, rarefaction waves will be reflected when the precompression 
wave of a subsonic flame reaches an area divergence, thus causing 
flame deceleration. However, experimental observations have shown 
that the transition from deflagration to detonation is enhanced rather 
than delayed as predicted on the basis of the simplified inviscid wave 
theory. 

From the foregoing physical description of the interaction phenom- 
ena with area changes (including obstructions) in an exothermic flow 
field, it can be foreseen that the cumulative effects will result in 
very complex flow patterns. On the basis of the wave -interaction 
analysis, the quantitative behavior of the interactions between shocks 
and area changes and between pressure waves and a flame discontinuity 
is sufficiently understood. The former has been investigated exten- 
sively by Oppenheim and Urtiev (Ref. 10), Rudinger (Ref. ll), Whitman 
(Ref. 12), Friedman (Ref. 13 ) and Rosciszewski (Ref* 14) and the latter 
by Rudinger (Ref. 15), and Chu (Ref. 1 6). In the present analysis, 
only the gas -dynamical aspect of the interaction between a flame and 
an area disturbance will be treated theoretically. A perturbation 
treatment following Chu f s approach (Refs. 1 6 and 17) will be employed. 
While a full scope of the subject matter cannot be explored by means 
of the perturbation theory, it nevertheless allows an insight to be 
obtained of the essential features of the interaction phenomena 
between a flame front and area disturbances. Furthermore, the 
complicated calculation procedure can be circumvented with this 
approach . 
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Symbols 


A 

AA 

a 

C P 

D 

a 

K 

L 


cross-sectional area of the duct 
finite change of the area 
speed of sound 

specific heat at constant pressure 

drag force parameter due to an obstruction 

pipe diameter 

correlation function 

detonation induction distance 



Mach number of the flame relative to unburned gas 


P 

Q 


E 

Vf d 

Re.. = — 
1 V 


S 

T 

T c 

u 


W = + Vf 


pressure 

heat release per unit mass of medium 
specific gas constant 

Reynolds number based on diameter and burning velocity 
entropy 

absolute temperature 
flame temperature 
gas velocity 
burning velocity 
flame speed 


5 


operator signifying infinitesimal change 
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p 


wv. 


pressure ratio 
temperature ratio 
ratio of specific heats 

ratio of specific gas constants across a flame front 
dens ity 


Subscripts 


1 conditions ahead of the flame 

2 conditions behind the flame 


B. Interaction Between a Subsonic Flame and a Cross-Sectional Area 
Change 

When two ducts of different cross-sectional areas are joined 
together, the flow field in which a flame propagates may be described 
approximately by a one -dimensional model containing a discontinuous 
change in cross section. This simplifying assumption presupposes the 
existence of a plane flame front which possesses an effective flame 
speed determined by the chemical kinetics and transport properties of 
the given mixture. Upon reaching the area discontinuity, the effective 
flame speed is influenced further by the interaction between the flame 
and the partially reflected precompression wave from the area discon- 
tinuity. In this analysis, attention is directed toward the gas- 
dynamical aspect of flame acceleration and deflagration-to-detonation 
transition during the initial build-up period for a flow of combus- 
tible gases in a duct with an area change. Thus, the flame front is 
assumed to approach the area discontinuity with a known "effective" 
flame speed at a low subsonic Mach number. 

Following the approach by Chu (Refs. l6 and 17) > a quasi -steady, 
one -dimensional flow model may be selected and the coordinate system 
may be oriented with respect to the flame front (Fig. 52). At the 
instant when the flame reaches the area discontinuity, infinitesimal 
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changes in properties , represented by 5, are assumed to take place 
both upstream and downstream of the flame front. As will be seen 
later, many higher-order terms of the perturbation quantities, such 


&p- bp 1 §p r 


as 


5W 

(but not necessarily can be neglected 


*1 1 2 

since they are small compared to unity. 

The conservation of mass across a flame front located at the 
area transition under the steady-state condition can be expressed as 


p i ( u i + v f) A = p 2 + v f^ d + &A ) d) 

However, it is assumed that small changes in properties occur as soon 

as the flame interacts with the area disturbance. In other words, 

o . d . u , u and v and so forth are now changed to p + 5p , 

K l 7 K 2 7 l 7 2 f _LX 

p 2 + bp 2’ ^2 + 5u 2> ^ V f + 5v f' etc * Tne continuit y equation (l) 
then becomes 

( p l + 5 p 1 ) (u x + 8u x + v f + 8v f ) A = 

( 2 ) 

( p 2 + 6 p 2 ) (u 2 + 5 u 2 + v f + 5v f ) (A + AA) 

In order to reduce equation (2) to a more convenient form, the 
following supplementary relations are introduced: 


u^ + v^ = W 


( 3 ) 


Su^ + 8v^ = 5W 


(L) 


Su^ + 5v f = 8W + (8u^ - Su^) 


( 5 ) 
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Ug + Vj. (u.t + Vj.) = — — ■ — i W 

^ t P 2 A + AA 1 f ^ aa 

A 


_ P 1 T p \ Pp (T 

where P = — , \ = — , (T= ~± and J- = JL. 

P 2 T 1 *2 Pi *€ 


( 6 ) 

(7) 


It is to be noted that the specific gas constant R as well as the 
specific heat c p and the ratio of specific heats 7, which are 

contained in the energy and entropy equations, generally have 
different values in the fluid media upstream and downstream of the 
flame front. However, their variations are regarded as being of 
little practical interest. 


With the 


imation f 1 + 


(■ 



the relations (3) through (7) and the approx- 

6p 2 

1 - n , the first-order expression for equation 


(2) becomes 


5u 2 - 5u-j_ 
W 



( 





+ 




(8) 


For a low subsonic flame for which the dynamic contribution is 
neglected, the momentum equation can be obtained by integrating Euler 1 s 
equation 


a lim =, 0 

dx 


(9) 
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which results in 


SPl Sp 2 P;L 

= with H = — = 1 

P P ^ P 

1 2 2 


( 10 ) 


Equation (10) implies that, for a small change in cross-sectional 
area, the pressure pulses ahead and behind a flame front are of 
the same type, i.e., either both compressive or both expansive. 
Furthermore , this equation has the same form as that for a constant- 
area duct. 

In a similar manner, the energy perturbation equation, ignoring 
the kinetic energy terms, can be expressed as follows: 

(in 


5Q = c 5T 2 - c BT-l 

*2 n 


Other conservation 
Entropy equation 


8p x 1 



laws in terms of the perturbation variables 


5p-j_ SS^ 



are : 


( 12 ) 


Equations of state 


5 P]l S P;l 



Sp 2 Sp 2 



BT-j 

h 

BT. 


(13) 

(14) 


Eliminating the temperature terms from the foregoing equations, the 
density variation can be expressed as follows : 
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Sp i 5C 2 


BP, 


1 2 


5Q , 6S 


V ‘Ri) 


1 


C P 

P 2 1 *1 


i . _2 ( 1 . c J!±l\ 

x c p n \ c p ? x ) 


Substitution of equation (15) into equation (8) yields 


6u 2 - 5u l / ^ AA\ 8W \ / AA\ 

w \ 1 + T/~T <jr " ( + Tj 


BP, 


(T 7 1 \ % V\ H / 


( 15 ) 


(16) 


BS, 


hfi-Oi iVi ♦ »\ 

<T\ C P, > J\ * ) \ T ! |<r\ »/ 


It is seen that for a constant-area duct, i.e., A A = 0, equation (l6) 
reduces to the expression obtained by Chu (Ref. 17). 

In order to analyze the type of pressure waves generated from 
the flame front located at the area disturbance, the perturbation 
quantities &W, 6Q and 5Sp in equation (l6) will have to be known 
or assumed. As a simplification, the effects of these variables on 
the generation of pressure waves are assumed to be additive, i.e.. 
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the overall effect can he treated as the sun of the individual 
effects. The influence of the change in flane speed on the flow 
characteristics under the condition of area variation will be 
singled out for detailed study. 

It is evident that additional conditions relating the pressure 
waves to the velocity perturbation and the flane speed to the area 
change are necessary for the analysis. The fomer can be fulfilled 
by introducing the characteristic relations whereas the latter can 
be supplenented partially on the basis of physical considerations. 
Since the flane front is located at the area transition from which 
pressure waves are generated, the characteristic equations for 
plane waves can be enployed to describe the pressure -velocity 
relationship in the two different constant -area ducts: 


bp-, Su-i 

— i + 7, — i = 0 


®1 


(17) 


&Po 5u o 

_±. - 7 P — = 0 


(18) 


Substitution of equations (17) and (l8) into equation (l6) while 
setting bQ = bS^ = 0 leads to 


bp-, 5p p \ 5W 

— = — = 1 7n M, — 

p p ) 1 1 w 

1 2 


f 


■ (‘ ' *)' 


( 19 ) 


HXfH' 


f i y 
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where M, = — is the Mach number of the flame front. It is seen 

-L a x 

that the denominator of equation (19) is always positive; therefore, 
whether the generated pressure waves are compressive or expansive 
will depend on the sign of the numerator. 

According to the present coordinate orientation (Fig. 52), a 
sudden area convergence is represented by AA > 0 and an area 
divergence is represented by A A < 0. For low subsonic flames, 

SW is positive for AA > 0 since the flow velocity increases in 
passing through an area constriction. Similarly, 5W is negative 
for AA < 0. Without any loss of generality, a simplified form of 
equation (19) can be obtained by assuming 7 = 7^ = 7^, = 1, and 

ignoring the second term in the denominator because of the low Mach 
number and the small area 


BPp Sp 2 7 Mjl 


riation: 


- 1 


1 + 


AA 

A 


6W 

vT 


( 20 ) 


Although in the case where AA > 0, the type of pressure waves depends 

on the relative magnitude between X and 1 + y , it appears that for 

most practical cases the bracketed term is positive for a small change 
in area. Other contributing factors which influence the determination 
of the type of pressure waves generated from the flame front are the 
changes in heat release and entropy when the area is varied. Mathe- 
matically, these effects can be analyzed by means of equation (l6) 
and the results are similar to those for a duct having a uniform 
cross section (Refs. 1 6 and 17). In general, the following conclusions 
hold: 


(i) For area convergence, AA > 0, 

(ii) For area divergence, AA < 0, 


Sp 


> 0 Compression 


^ < 0 Expansion 
P 


( 21 ) 
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The effect of the generated pressure waves on further flame 
characteristics can be predicted from physical considerations. If 
compression waves are produced, the rearward propagating wave will 
be reflected from the closed end of the duct and this wave will 
catch up with the flame front thereby creating a new set of trans- 
mitted and reflected compression waves. Thus, the flame will be 
accelerated further and intermittent pressure-wave reinforcement 
will take place leading to the eventual formation of strong shock 
waves capable of auto -ignition. On the other hand, expansion waves 
will produce opposite effects which tend to attenuate the flame and 
shock intensities. 

While the influence of area variation on the formation of 
detonation waves can be analyzed qualitatively based on the pertur- 
bation theory, this simplified fluid model selected herein may not 
be adequate, under certain circumstances, to explain the complex 
phenomena that actually exist. This aspect will be discussed 
further in a later section. 


C. Interaction Between a Subsonic Flame and a Uniform Rod-Insert 
in a Constant-Area Duct 


At the outset, this problem appears to be more complicated 
than that involving a single area discontinuity considered previously 
since the presence of a rod-insert not only modifies the cross- 
sectional area but also leads to skin friction on the surface of the 
rod. However, the perturbation theory still can be applied to a 
limited extent if the obstruction (the rod-insert) is assumed to be 
represented by a constant skin-friction parameter. The problem 
then reduces to one with an area constriction as well as that 
involving force and energy modifications (Fig. 52). The continuity 
equation remains the same as equation (2) or (8) except that the 


pressure 

friction 
based on 


ratio is of the form £ = 



Pi 

P 1 + D 


where D is a skin- 


parameter due to the obstruction. The momentum equation 
the one -dimens ional assumption then can be expressed as 


P 1 + 6 P;l + D = p 2 + Sp 2 


( 22 ) 
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or 
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5p 2 

Therefore, the reflected wave represented by - — is of the same type 

Bp, 2 

as the transmitted wave , — — , but the reflected wave is weaker in 

P 1 


strength by a factor € . The energy perturbation equation still can 
be represented by equation (11) provided that the heat addition term, 
q, includes the effects due to chemical reaction as well as those 
resulting from viscous dissipation. 


Proceding with the analysis as before, the condition which the 
perturbed flame must satisfy is found to be: 



On the basis of the argument that 6W increases for an area 
decrease due to the presence of the rod, the simplified expression 
relating the generated pressure waves and the change in flame 
speed ^becomes 
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5pi 


P 


1 


1 ^P2 

€ p 2 


? M i 

j*7 + i 




1 + 


AA 

A 



(25) 


Equation (25) shows that if £ is not too small and if the diameter of 
the rod is not too large, compression waves will he generated when 
the flame front reaches the rod-insert. 

The limitation of the perturbation treatment must be recognized 
since far-reaching simplifications have been made in the analysis. 

The validity of a one -dimensional fluid model is questionable if 
the rod-insert has a relatively large diameter. Under this circum- 
stance, the boundary -layer effect is expected to be quite pronounced 
and the influence of turbulence generation cannot be disregarded 
totally. 


D. Empirical Relations 

In previous studies (see Ref. 6), the transition phenomena from 
deflagration to detonation have been analyzed empirically. Of 
particular interest was the measurement of the detonation induction 
distance, a length parameter characterizing the development of deto- 
nation waves. It has been found that a correlation function defined 
as 


K = Re-L 



(26) 


exhibits an almost unique relationship with the statistically 
determined detonation induction distance for a given constant-area 
duct. The empirical results also show that K is an almost universal 
function in the sense that the functional relationship holds for a 
number of non -hydrocarbon systems regardless of the pressure and 
fuel concentration. Because of • insufficient data, this correlation 
function, K, has not been tested previously for ducts of different 
diameters. Recently, however, a substantial number of measurements 
of the detonation induction distance for various hydrogen-oxygen 
mixtures in constant-area ducts of diameters of 5 cm and 7*9 cm have 
been completed. It has been found that a dimensionless length, 
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— , vhere L is the induction distance, is a more suitable parameter 


for empirical studies. 

1/4 


between 


(i) 


1/4 


Figure 53 depicts the functional relationship 
for various hydro gen -oxygen mixtures at 


different pressures with d = 5 cm, 7*9 cm, and also the results 
obtained previously for hydrogen -oxy gen and carbon monoxide - oxygen 
mixtures (Ref. 2) for a pipe with d = 1.5 cm. The numerical eval- 
uations are given in Tables 3 through 6. Although the results show- 
some scattering on the graph, the deviation of data from the average 
curve is not serious enough to discount the universal property of the 
correlation function K. One possible cause for the scattering of 
data comes from the calculation of K for which some of the values of 
the burning velocity, v^., were extrapolated from experimental curves, 

and the combustion temperature, T c , was based on theoretical calcu- 
lations rather than from an experimental determination (Ref. 6). The 
advantage of this approach lies in the fact that the formation of 
detonation waves for a non-hydrocarbon mixture in a constant-area 
pipe can be predicted, within some degree of accuracy, by the 
correlation function K which depends only on the properties of the 
initially subsonic flame. 


The empirical analysis discussed above has been extended to 
problems involving a non-uniform cross-section a l area and the presence 
of a rod-insert in a constant-area pipe. A method of superposition, 
similar to the determination of the wake function in turbulent flow 
theory, has been attempted in order to account for the influence of 
an area variation or an obstruction on the transition mechanism. 
However, in view of some of the peculiarities of the experimental 
results (see Table l), a suitable correlation function has not yet 
been found. 


E. Discussion 

On the basis of the one -dimensional perturbation theory, it has 
been shown that the interaction between a flame front and an area 
change results in the generation of pressure waves propagating both 
upstream and downstream of the flame front. In general, compression 
waves are generated for area convergence and expansion waves are 
generated for area divergence. If the disturbance is a rod-insert, 
compression waves will be produced. 

From equation (l6) it is seen that the pressure perturbation 
depends on the changes in flame speed, heating value, and entropy 
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under the influence of an area discontinuity • As evidenced from the 
present mathematical procedure, the effect of an area change on wave 
generation cannot be predicted directly. However, it can be conceived 
that the non -uniformity of the geometric configuration of a duct 
introduces an additional condition which affects the variation of the 
physical quantities, such as the flame propagation rate. The latter, 
in turn, will determine the type of the generated pressure waves. For 
example, upon reaching an area convergence, the unburned gas ahead of 
the flame front will accelerate and the reaction zone will be thick- 
ened. Therefore, an increase in flame propagation rate will result 
and 5W will be positive, thereby leading to the generation of compres- 
sion waves as discussed previously. Mathematically speaking, the 
present problem can be regarded as that for a constant-area duct 
for which the signs of SW, 5Q, and so forth, have been predetermined 
by the area effects. It should be mentioned that the change of flame 
speed through a duct with non-uniform cross section has been predicted 
purely on physical grounds since a rigorous mathematical analysis of 
this problem is not yet available. 

The actual flow field for a combustion wave propagating in a 
duct with an area discontinuity is a very complex one. The present 
analysis is concerned only with the interaction between the flame and 
the area discontinuity under low subsonic conditions. However, as 
soon as combustion is initiated from the closed end, a precompression 
wave will be present in order to satisfy the boundary conditions . 

This compression wave will reach the area discontinuity before the 
flame front and this wave will be partially transmitted and partially 
reflected. With the assumption that the unburned gas flow is subsonic, 
the reflected wave will propagate rearward instead of being "swallowed" 
by the area convergence or divergence as that under high supersonic 
conditions . 

For area convergence, the reflected wave is compressive (weak 
shock) which will interact with the ensuing flame front. As discussed 
by Rudinger (Ref. 15), the head-on collision between a shock and a 
flame front will distort the flame surface and the indentation will 
develop quickly into a spike. The flame area increases greatly and 
strong turbulence results. Therefore, flame acceleration occurs which 
will create further compression waves from the flame front. Moreover, 
as a result of the shock-flame interaction, rearward moving trans- 
mitted shock waves will be reflected from the closed end and will 
overtake the flame causing further modification of the flame front. 

As the flame enters into the region of a different cross-sectional 
area, the interplay among the area disturbance, the flame, and the 
shocks will result in additional waves being generated. In spite of 
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the complex flow field that arises, it can be foreseen that these 
effects will lead to the reinforcement of the intensity of the 
precompression wave and to a continuing acceleration of the flame. 

The transition from deflagration to detonation therefore is expected 
to be enhanced. 

For area divergence, it appears that because the precompression 
wave is reflected from the area disturbance as rarefaction waves, 
the intensity of the flame may be weakened. However, for a sudden 
area expansion, vortices are generated and flow separation will 
result. Thus, the mechanism of turbulence generation will compen- 
sate for the weakening of the shock (precompression wave) intensity 
due to the wave -interaction phenomena. The final outcome still is 
the enhancement of the transition process from deflagration to 
detonation. Indeed, the experimental measurement of the detonation 
induction distance for hydrogen -oxygen combustion waves through a 
duct with sudden area enlargement clearly confirms this conclusion. 

If a uniform rod is inserted symmetrically in a constant -area 
duct, the effects are manifold: 

(1) The cross-sectional area is reduced as if two uniform ducts 
of different cross-sectional area were joined together, 

(2) The precompression wave as well as the flame will impinge 
on the obstruction (the front surface of the rod-insert) 
causing a series of wave -interaction phenomena, 

(3) The rod-insert furnishes an additional source of real gas 
effects, such as the effects from viscous stresses and heat 
transfer . 

From a qualitative viewpoint, the transition phenomena can be 
analyzed or at least predicted by considering the individual effects 
listed above. The wave -interaction aspect of the problem follows 
the same procedure as before, but the generation of turbulence is 
determined partially by the shape and the diameter of the rod-insert. 
For larger diameter rod-inserts, the heat-sink effect becomes more 
significant, and the one -dimensional flow model may be inadequate 
to describe the flow field fully. 

Whether the disturbance is caused by area discontinuity or by a 
rod-insert, it can be concluded that the transition from deflagration 
to detonation generally will be enhanced as compared with that for a 
constant-area duct. Without sufficient quantitative information about 
the effects of area change and obstruction on flow disturbances, the 
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problem is not amenable to empirical analysis at the present time. 

Also some difficulties were encountered in the determination of the 
detonation induction distances for a larger diameter rod-insert 
(d = 4.76 cm) for which two sets of additional supports had to be 
used in order to center the rod. The measured velocities near these 
supports were found to have large values locally which indicate a 
different experimental condition from that without these supports . 

The task of determining the induction distance accurately, therefore, 
was quite difficult. 

Although some understanding can be gained by employing the method 
of characteristics to analyze the interaction phenomena between a 
flame front and an area disturbance or obstruction, the calculated 
flow patterns still are very likely to deviate from the real ones. 

In particular, the important role of the turbulence generation under 
these circumstances cannot be incorporated easily to the graphical 
and iteration procedures. 


CONCLUSIONS 


When a combustion wave reaches an area discontinuity or a rod- 
insert, pressure waves will be generated. These generated pressure 
waves will induce subsequent wave -interaction phenomena which will 
strongly influence the transition process from deflagration to 
detonation. On the basis of the small perturbation analysis, it is 
shown that the effects of the area disturbances generally will 
enhance the transition process, as verified by measurements of the 
detonation induction distances . 

The results of the measurements of detonation induction distances 
with various rod-inserts and geometrical configurations are included 
in this analysis. It is evident that the mechanism of turbulence 
generation is quite important when these disturbances sire present 
in the flow field. Without consideration of the effects of turbulence, 
the transition mechanism for a flame propagating in a duct with non- 
uniform cross section cannot be explained satisfactorily. 

The simplified flow model selected for theoretical analysis of 
the interaction phenomena between a low subsonic combustion wave and 
area disturbances permits a qualitative understanding of the complicated 
problem. Although a general trend of the transition process from 
def lag ration to detonation due to area effects (including an obstruc- 
tion) can be predicted according to the perturbation theory, a great 
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improvement in the analysis can be achieved if viscous effects are 
incorporated therein. 

Some of the data from the measurements of the detonation induction 
distances (Table l) show surprising results. For example, the induction 
distance for a stoichiometric hydrogen-oxygen mixture at 5 atmospheres 
initial pressure varies in a somewhat inconsistent manner with the 
change in the diameter of the rod-insert. At 5 atmospheres pressure, 
the induction distances increase with the amount of blockage until the 
rod-insert with 60 per cent blockage is employed whereupon the induction 
distance decreases drastically. Because of the limited amount of data 
available and until a satisfactory explanation of this peculiarity is 
given, it is difficult to develop an empirical formula similar to that 
for a constant-area duct. 


The Ohio State University 

Rocket Research Laboratory 

Columbus 10, Ohio August 1, 1962 
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TABLE 2. COMPARISON OF DETONATION INDUCTION DISTANCES FOR HYDROGEN-OXYGEN 
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TABLE 4. DETONATION INDUCTION DISTANCES FOR THE RIGHT ANGLE 
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TABLE 5* 


EVALUATION OF 



1 

4 


FOR 


VARIOUS HYDROGEN-OXYGEN MIXTURES (d = 1.5 cm) 


Mol 

Per Cent 
Hydrogen 

Pi 

(atm) 

L 

(cm) 

K 

/L \ 
\ d / 

(if 

(f 

40.0 

1 

102 

276.9 

68 

2.87 

4.07 


5 

41 

3076.0 

27.3 

2.28 

7.42 


10 

36 

9409.5 

24 

2.22 

9.85 


25 


44315.1 

“ 

- 

14.5 

50.0 

l 

66 

481.0 

44 

2.58 

4.68 


5 

20 

5508.4 } 

13.3 

1.91 

8.61 


10 

13.5 

16714.6 

9 

1.73 

11.37 


25 

5 

79839.1 

3.4 

1.36 

16.79 

6o.o 

1 

6l 

1026.8 

40.6 

2.52 

5.66 


5 

16 

9340.2 

10.67 

1.81 

9.82 


10 

13 

25663.6 

8 .68 

1.72 

12.65 


25 

5 

108187.6 

3.4 

1.36 

18.13 

66.7 

1 

_ 

1278.2 


- 

5.97 


5 

19 

9565 .8 

12.66 

I.89 

9.88 


10 

12 

24767.0 

8.0 

1.68 

12.54 


25 

2.5 

95239.3 

1.67 

1.14 

17.58 

70.0 

1 

mm 

1222.4 

60.7 

2.79 

5.91 


5 

m 

7510.4 

13.3 

1.91 

9-3 


10 


18317.5 

8.68 

I.72 

11.64 


25 

iiiti 

69281.4 

2.67 

1.28 

16.23 

8o.o 

1 

156 

290.6 

104 


4.18 


5 

53 

1743.8 

35-35 


6.45 


10 

45 

3869.5 

30 


7.9 


25 

31 

12113.0 

20.65 


10.52 


Values of L and K are obtained from Ref. 6 
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TABLE 6. EVALUATION Of|k^ AND ^ | j ^ FOR 


VARIOUS HYDROGEN-OXYGEN MIXTURES (d = 5.0 cm) 


Mol 

Per Cent 
Hydrogen 

P 1 

(atm) 

L 

(cm) 

K 

■ 


H 1 

45.0 

1 

100 

1315 

20 

2.11 

6.02 

60.0 

1 

50 

3420 

10 

1.78 

7.65 

66. 7 

1 

70 

4260 

i4 

1.93 

8.07 

70.0 

1 

70 

4o8o 

l4 

1.93 

7.96 

75-0 

1 

166 

2135 

33.2 

2.4 

6.80 

80.0 

1 

290 

987 

58 

2.76 

1 

5.60 


Values of L are obtained from Ref. 2. 
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TABLE 7 . 


EVALUATION OF 



1 

4 


FOR 


VARIOUS HYDROGEN -OXYGEN MIXTURES (d = 7.9 cm) 


Mol 

Per Cent 
Hydrogen 

45.0 

66.7 

75-0 


Pi 

(atm) 

L 

(cm) 

K 

Q 


■ 

1 

186 

1922 

23.55 

2.2 

6.62 

5 

86 

21900 

10.9 

1.82 

12.16 

1 

137 

6720 

17.33 

2.04 

9.05 

5 

37 

50300 

4.7 

1.47 

14.97 

1 

207 

3365 

26.2 

2.26 

7.61 

5 

90 

24200 

11.4 

1.84 

12.48 
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TABLE 8. 


EVALUATION OF 



1 

4 


FOR 


VARIOUS CARBON MONOXIDE - OXYGEN MIXTURES (d = 1.5 cm) 


Mol 

Per Cent 
Carbon Monoxide 

Pi 

(atm) 

L 

(cm) 

K 

■ 

■ 

■ 

45. 0 

5 

156 

57.8 

104 

3-20 

2.76 


10 

l4l 

128.7 

9^.0 

3-12 

3-37 


25 

107 

388.6 

71.3 

2.91 

4.1*4 

60.0 

5 

129 

123.2 

86.0 

3.04 

3.33 


10 

122 

318.2 

81.3 

3.00 

4.26 


25 

74 

1294.5 

49.3 

2.65 

6.00 

66.7 

5 

125 

152.7 

83.3 

3.02 

3.51 


10 

117 

392.4 

79.O 

2.98 

4.46 


25 

6l 

1703.1 

Uo. 7 

2.52 

6.42 

75.0 

5 

143 

172.3 

95-3 

3.12 

3.62 


10 

125 

1*1.8 

83.3 

2.65 

4.58 


25 

114 

1917.8 

76.0 

2.95 

6.61 

80.0 

5 

154 

147.3 

102.7 

3.18 

3.48 


10 

145 

356.3 

96.7 

3.14 

4.35 


25 

126 

1349.7 

81*. 0 

3.02 

6.05 


Values of L and K are obtained from Ref. 6. 
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FIG. 1 PHOTOGRAPH OF DETONATION TUBE 






FIG. 3 PHOTOGRAPH OF FLOW CONTROL SYSTEM 
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FIG. 5 PHOTOGRAPH OF CHRONOOIAPH SYSTEM 





FIG, 6 DESCRIPTIVE VIEW OF CYLINDRICAL INSERT, DIAMETER = 0.635 CM 
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FIG. 7 DESCRIPTIVE VIEW OF CYLINDRICAL INSERT, DIAMETER = I.588 CM 
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DESCRIPTIVE VIEW OF CYLINDRICAL INSERT, DIAMETER = 4.763 CM 





FIG. 9 DESCRIPTIVE VIEW OF DIVERGENT STEPPED-WALL INSERT 
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FIG. 11 RATE OF FLAME PROPAGATION IN Hg - 0 2 MIXTURE; MOLE # H 
V± = 1 ATM, t ± = hO°C, CLEAR TUBE. 
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FIG. 12 RATE OF FLAME PROPAGATION IN Hg - Og MIXTURE; MOLE % Hg 
p, = 1 ATM, ti = 40°C, CLEAR TUBE. 
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FIG. 13 RATE OF FLAME PROPAGATION IN Hg - Og MIXTURE; MOLE j> Hg = 75 
p. = 1 ATM, t n . = 40°C, CLEAR TUBE. 
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14 RATE OF FLAME PROPAGATION IN Hg - Og MIXTURE; MOLE £ Hg 
p., = 5 ATM, t t = 40°C, CLEAR TUBE. 
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FIG. 1 6 RATE OF FLAME PROPAGATION IN Hg - 0 2 MIXTURE; MOLE j> Hg 
p. = 5 AIM, t, = 40°C, CLEAR TUBE. 
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FIG. 17 RATE OF FLAME PROPAGATION IN H 2 - 0 2 MIXTURE; MOLE $ Hg = 66 . 67 , 
Pj_ = 1 ATM, t ± = 40°C, 0.635 -CM DIAMETER INSERT. 
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FIG. 18 RATE OF FLAME PROPAGATION IN Hg - 0 2 MIXTURE; MOLE $ Hg = 66 . 67 . 
p, = 5 a™, t -5 = 40°C, O. 635 -CM DIAMETER INSERT. 
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FIG. 19 BATE OF FLAME PROPAGATION IN Hg - 0 2 MIXTURE; MOLE $ Kg 
p. = 1 ATM, = 40°C, 1.588 -CM DIAMETER INSERT. 
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FIG. 20 RATE OF FLAME PROPAGATION IN Hg - 0 2 MIXTURE; MOLE $ Hg 
p. = 5 ATM, tp = 40°C, I. 588 -CM DIAMETER INSERT. 
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FIG. 22 RATE OF FLAME PROPAGATION IN Hg - Og MIXTURE; MOLE % EL = 66.67, 
P ± = 1 AIM, t i = 40°C, 4 . 763 -CM DIAMETER INSERT. 
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24 RATE OF FLAME PROPAGATION IN Hg - 0 2 MIXTURE; MOLE <f> Hg = 45 
p. = 5 ATM, t, = 40°C, 4. 763 -CM DIAMETER INSERT. 
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FIG. 25 RATE OF FLAME PROPAGATION IN Hg - 0 2 MIXTURE; MOLE <f> IL> = 66.67 
P i = 5 ATM, t ± = kO°C, 4.763-CM DIAMETER INSERT. 



FIG. 26 RATE OF FLAME PROPAGATION IN Hg - 0 £ MIXTURE; MOLE £ Hp = 75, 
p i *= 5 AIM, = 40°C, 4.763-CM DIAMETER INSERT. 
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FIG. 27 DETONATION INDUCTION DISTANCES OF HYDROGEN -OXYGEN MIXTURES 
BASED ON MAXIMUM FLAME PROPAGATION RATES FOR THE CLEAR TUBE 
AND 0.635, 1.588 AND 4.763-CM DIAMETER INSERTS. 
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FIG. 28 DETONATION INDUCTION DISTANCES OF HYDROGEN-OXYGEN MIXTURES 
BASED ON AVERAGE FLAME PROPAGATION RATES FOR THE CLEAR TUBE 
AND 0.635, 1.588 AND 4.763-CM DIAMETER INSERTS. 
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29 RATE OF FLAME PROPAGATION IN Hq - Og MIXTURE; MOLE i> Hg 
p. = 1 AIM, t.; = 40°C, DIVERGENT STEPPED-WALL INSERT. 
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30 RATE OF FLAME PROPAGATION IN Hg - Og MIXTURE; MOLE # Hg 
Pj = 1 ATM, -4 = 40°C, DIVERGENT STEPPED-WALL INSERT. 
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FIG. 31 RATE OF FLAME PROPAGATION IN Hg - Og MIXTURE; MOLE % Hg 
p, = 1 ATM, t,- = hO°C, DIVERGENT STEPPED -WALL INSERT. 
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FIG. 32 DETONATION INDUCTION DISTANCES OF Hg - 0 2 MIXTURES BASED ON 
MAXIMUM FLAME PROPAGATION RATES WITH A DIVERGENT STEPPED - 
WALL INSERT. 
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FIG. 33 DETONATION INDUCTION DISTANCES OF Hg - 0 2 MIXTURES BASED ON 

AVERAGE FLAME PROPAGATION RATES WITH A DIVERGENT STEPPED - 
WALL INSERT. 
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FIG. 35 RATE OF FLAME PROPAGATION IN FL, - 0 2 MIXTURE; MOLE j> Ho = 66.67 
p = 1 AIT-!, t i = 40°C, 90° -BEND 66 CM FROM IGNITOP. 
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FIG. 36 RATE OF FLAME PROPAGATION IN Hg - 0 2 MIXTURE; MOLE $ Hp = 75, 
P. = 1 ATM, t. = 40°C, 90° -BEND 66 CM FROM IGNITOR. 
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FIG. 38 RATE OF FLAME PROPAGATION IN Hg - Og MIXTURE; MOLE 3 Hg = 80 . 0 'i 
p. = 5 ATM, t, = 40°C , 90° -BEND 66 CM FROM IGNITOR. 
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FIG. 39 RATE OF FLAME PROPAGATION IN H2 - 0 2 MIXTURE; MOLE # IL =75, 
p j[ = 5 ATM, t i = 40°C , 90° -BEND 66 CM FROM IGNITOR. 



FIG. 40 RATE OF FLAME PROPAGATION IN Hg - 0 2 MIXTURE; MOLE Hg = 
P i = 10 ATM, t i = 40°C, 90° -BEND 66 CM FROM IGNITOR. 
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FIG. 4l RATE OF FLAME PROPAGATION IN Hg - Og MIXTURE; MOLE $ Hg = 66.67 
p i = 10 ATM, t ± = 40°C, 90° -BEND 66 CM FROM IGNITOR. 
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FIG. 42 RATE OF FLAME PROPAGATION IN Hg - 0 2 MIXTURE; MOLE # Hg =75, 
p. = 10 ATM, t, = 40°C, 90° -BEND 66 CM FROM IGNITOR. 
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FIG. 43 RATE OF FLAME PROPAGATION IN Hg - 0 2 MIXTURE; MOLE 
p, = 1 ATM, = 40°C, 90 ° -BEND 15 CM FROM IGNITOR 
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FIG. 44 RATE OF FLAME PROPAGATION IN Hg - 0 2 MIXTURE; MOLE j> = 66.67 
p, = 1 AIM, t,- = 40°C, 90° -BEND 15 CM FROM IGNITOR. 
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FIG. k5 RATE OF FLAME PROPAGATION IN Hg - 0 2 MIXTURE; MOLE $ Hg 
p = 1 ATM, t i = 40°C, 90° -BEND 15 CM FROM IGNITOR. 
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FIG. bj RATE OF FLAME PROPAGATION IN Hg - Og MIXTURE; MOLE $ Hg = 66.67 
p. = 5 ATM, t. = 40°C, 90° -BEND 15 CM FROM IGNITOR. 
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FIG. bQ RATE OF FLAME PROPAGATION IN Hg - Og MIXTURE; MOLE <f> Hg 
V, = 5 AM, t ± = 40°C , 90° -BEND 15 CM FROM IGNITOR. 
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RATE OF FLAME PROPAGATION IN H 2 - 0 ? MIXTURE; MOLE Hg 
p = 10 AIM, t i = 40°C, 90° -BEND 15 CM FROM IGNITOR. 
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50 RATE OF FLAME PROPAGATION IN H 2 - O 2 MIXTURE; MOLE <f> Hg 
p, = 10 ATM, t ± = LO°C, 90° -BEND 15 CM FROM IGNITOR. 
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FIG. 51 RATE OF FLAME PROPAGATION IN Hg - 0 2 MIXTURE; MOLE # Hg = 75, 
p i = 10 ATM, t ± = 40°C, 90° -BEND 15 CM FROM IGNITOR. 
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FIG. 52 COORDINATE SYSTEM ORIENTATION FOR A FLAME IN A DUCT 
WITH AN AREA DISCONTINUITY OR AN OBSTRUCTION. 
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FIG. 53 DIMENSIONLESS INDUCTION DISTANCE AS A FUNCTION OF THE 
CORRELATION FUNCTION K. 
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